arXiv:1005.2366vl  [cond-mat.quant-gas]  13  May  2010 


Competition  between  pairing  and  ferromagnetic  instabilities  in  ultracold  Fermi  gases  near 

Feshbach  resonances 


David  Pekker1,  Mehrtash  Babadi1,  Rajdeep  Sensarma2,  Nikolaj  Zinner1,3,  Lode  Pollet1,  Martin  W.  Zwierlein4,  Eugene  Demler1 
1  Physics  Department,  Harvard  University,  Cambridge,  Massachusetts  02138,  USA 
2  Condensed  Matter  Theory  Center,  University  of  Maryland,  College  Park,  Maryland  20742,  USA 
3  Department  of  Physics  and  Astronomy,  Aarhus  University,  Aarhus,  DK-8000,  Denmark 
4  MIT-Harvard  Center  for  Ultracolcl  Atoms,  Research  Laboratory  of  Electronics, 
and  Department  of  Physics,  Cambridge,  MA  02139,  USA 

We  study  the  quench  dynamics  of  a  two-component  ultracold  Fermi  gas  from  the  weak  into  the  strong  inter¬ 
action  regime,  where  the  short  time  dynamics  are  governed  by  the  exponential  growth  rate  of  unstable  collective 
modes.  We  obtain  an  effective  interaction  that  takes  into  account  both  Pauli  blocking  and  the  energy  depen¬ 
dence  of  the  scattering  amplitude  near  a  Feshbach  resonance.  Using  this  interaction  we  analyze  the  competing 
instabilities  towards  Stoner  ferromagnetism  and  pairing. 


Ferromagnetism  in  itinerant  Fermions  is  a  prime  example 
of  a  strongly  interacting  system.  Most  theoretical  treatments 
rely  on  a  mean-field  Stoner  criterion  m,  but  whether  this  ar¬ 
gument  applies  beyond  mean-field  remains  an  open  problem. 
It  is  known  that  the  existence  of  the  Stoner  instability  is  very 
sensitive  to  the  details  of  band  structure  and  interactions  ESI, 
however  how  to  account  for  these  details  in  realistic  systems 
remains  poorly  understood.  Exploring  the  Stoner  instability 
with  ultracold  atoms  has  recently  attracted  considerable  atten¬ 
tion.  Following  theoretical  proposals  a  ,  the  MIT  group  made 
use  of  the  tunability  {6|  and  slow  time  scales  lf7UT0l  of  ultra¬ 
cold  atom  systems  to  study  the  Stoner  instability  HD.  Sig¬ 
natures  compatible  with  ferromagnetism,  as  understood  from 
mean-field  theory  flZl.  were  observed  in  experiments:  a  max¬ 
imum  in  cloud  size,  a  minimum  in  kinetic  energy  and  a  maxi¬ 
mum  in  atomic  losses  at  the  transition.  However,  no  magnetic 
domains  were  resolved. 

An  important  aspect  of  the  MIT  experiments  is  that  they 
were  done  dynamically:  the  Fermi  gas  was  originally  prepared 
with  weak  interactions  and  then  the  interactions  were  ramped 
to  the  strongly  (repulsive)  regime.  Dynamic  rather  than  adi¬ 
abatic  preparation  was  used  in  order  to  avoid  production  of 
molecules.  This  raises  the  question  of  what  are  the  dominant 
instabilities  of  the  Fermi  gas  in  the  vicinity  of  a  Feshbach  res¬ 
onance. 

Naively,  one  would  expect  that  on  the  BEC-side,  molecule 
production  is  slow,  as  it  requires  a  three-body  process.  There¬ 
fore,  instability  towards  Stoner  ferromagnetism  would  dom¬ 
inate  over  the  instability  toward  molecule  production.  In 
this  picture,  one  would  expect  that  quenches  to  the  attrac¬ 
tive  (BCS)  regime  always  yield  an  instability  towards  pairing, 
whereas  quenches  to  the  repulsive  (BEC)  regime  an  instability 
towards  ferromagnetism  for  sufficiently  strong  interactions. 

In  this  Letter,  we  argue  that  this  picture,  which  was  used  to 
interpret  the  MIT  experiments,  is  incomplete.  Near  the  Fes¬ 
hbach  resonance,  even  on  the  BEC  side,  pair  production  re¬ 
mains  a  fast  two-body  process  as  long  as  the  Fermi  sea  can 
absorb  the  molecular  binding  energy.  As  a  result,  near  the 
Feshbach  resonance,  both  on  the  BEC  and  the  BCS  side,  the 
pairing  and  the  Stoner  instabilities  compete  directly.  We  now 
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FIG.  1 :  Growth  rate  of  the  pairing  and  Stoner  ferromagnetic  insta¬ 
bilities  after  a  quench  as  a  function  of  the  final  interaction  strength 
1  /kpa.  Final  interactions  with  negative  (positive)  values  of  1  /kpa 
correspond  to  the  BCS  (BEC)  side  of  the  Feshbach  resonance.  The 
Stoner  instability  simultaneously  occurs  in  multiple  channels.  The 
most  unstable  channel  is  indicated  by  the  solid  red  line,  the  others 
by  dashed  red  lines.  The  “RPA  Stoner”  instability  corresponds  to  the 
RPA  result  with  bare  as  opposed  to  Cooperon-mediated  interaction 
(see  text  and  Ref.  cm). 

Inset:  Schematic  diagram  of  the  pair  creation  process  showing  the 
binding  energy  (spring)  being  absorbed  by  the  Fermi  sea  (arrows). 


discuss  these  instabilities  and  their  competition  in  detail. 

We  start  by  describing  the  inter-atomic  interactions.  A  Fes¬ 
hbach  resonance  enables  tunable  interactions  between  ultra¬ 
cold  atoms  by  coupling  the  collision  partners  to  a  molecular 
state  with  different  magnetic  moment.  For  broad  resonances, 
where  the  coupling  is  much  larger  than  the  Fermi  energy,  this 
can  be  modelled  by  a  single  collision  channel  that  supports 
one  shallow  bound  state  m.  An  often  used,  but  patholog¬ 
ical  choice,  is  to  describe  repulsive  interactions  with  a  hard- 
sphere  pseudo-potential.  Instead,  one  should  use  the  full  T- 
matrix  that  includes  the  molecular  bound  state  ED.  Although 
at  low  energies  the  scattering  amplitude  from  the  hard-sphere 
potential  and  the  T-matrix  match,  at  higher  energies  compara¬ 
ble  to  the  molecular  binding  energy,  they  do  not.  Specifically, 
in  the  strong  interaction  regime  where  the  Stoner  instability 
is  expected  to  occur,  the  Fermi  energy  is  comparable  to  the 
binding  energy  of  a  molecule  in  vacuum,  precluding  the  use 
of  the  hard-sphere  potential. 
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In  light  of  this  remark,  we  study  the  initial  dynamics  of 
the  collective  modes  of  a  Fermionic  system  after  a  sudden 
quench,  taking  the  Cooperon  (full  T-matrix  and  Pauli  block¬ 
ing)  into  account.  We  focus  on  the  case  of  a  sudden  quench, 
as  it  is  simpler  and  captures  the  essential  physics  of  the  insta¬ 
bility  of  the  Fermi  surface.  Extensions  to  finite  rate  quenches 
are  discussed  in  Ref.  tm  Our  main  findings  are  summa¬ 
rized  in  Fig.[l]and  are:  (1)  We  find  that  with  the  full  T-matrix 
the  Stoner  instability  survives  with  a  finite  growth  rate  in  the 
range  —0.2  <  kFa  <  1,  where  a  is  the  scattering  length  and 
kp  is  the  Fermi  momentum.  In  contrast,  using  bare  interac¬ 
tions  ED  results  in  an  unphysical  divergence  of  the  growth 
rate  at  unitarity  and  no  magnetic  instabilities  on  the  BCS  side 
(see  Fig.  [TJ.  (2)  The  pairing  instability  persists  on  the  BEC 
side,  where  it  competes  with  the  Stoner  instability.  (3)  Within 
our  approximations,  the  pairing  instability  is  stronger. 

At  first  sight,  the  survival  of  the  pairing  and  Stoner  instabil¬ 
ities  on  the  wrong  side  of  the  resonance  is  quite  remarkable. 
However,  both  can  be  understood  by  taking  into  account  the 
presence  of  the  Fermi  sea.  On  the  BEC  side,  due  to  Pauli 
blocking,  the  binding  energy  of  the  pair-like  molecule  can  be 
absorbed  by  the  two  holes  that  are  left  behind  (see  the  inset 
of  Fig.  [TJ.  Thus,  the  two-body  pairing  process  becomes  for¬ 
bidden  when  the  binding  energy  ~  1/tocj2  exceeds  the  maxi¬ 
mum  energy  that  can  be  absorbed  by  the  holes  ~  kp/m(rn  is 
the  Fermion  mass,  a  the  scattering  length,  fcp  the  Fermi  mo¬ 
mentum,  and  throughout  this  Letter  we  use  the  units  in  which 
H  =  1).  On  the  BCS  side,  although  interactions  at  low  ener¬ 
gies  are  indeed  attractive,  the  same  is  not  true  at  high  energies. 
As  the  Stoner  instability  involves  all  scattering  energies  up  to 
the  Fermi  energy,  it  is  natural  that  it  can  persist  around  unitar¬ 
ity,  even  on  the  BCS  side. 

Formalism  -  We  consider  a  system  of  interacting  Fermions 
described  by  the  Hamiltonian: 

H  =  X^cLcka  +  J  ddicU{t.r (1) 

k,cr 

where  cj.( ca)  are  the  Fermion  creation  (annihilation)  opera¬ 
tors  with  spin  er,  £k(T  =  k2 /2m  —  pa,  are  the  chemical 
potentials,  and  U(t,r  —  r')  is  the  time  dependent  pseudo¬ 
potential  that  describes  the  inter-atomic  interaction.  We  fo¬ 
cus  on  the  instantaneous  quench  limit,  in  which  the  coupling 
U  changes  from  a  negligible  initial  value  Ul  to  a  final  value 
Uf  at  time  t  =  0.  In  this  limit,  we  can  describe  short  time 
dynamics  of  a  collective  mode  at  momentum  q  using  the  cor¬ 
responding  susceptibility,  xq(wq;  U /),  evaluated  with  final  in¬ 
teractions  but  initial  Fermionic  configuration  ED  Ei  In  par¬ 
ticular,  if  Xq(wq!  Uf)  has  poles  at  tcq  =  Hq  +  iAq  in  the 
upper  half  of  the  complex  plane,  then  fluctuations  that  occur 
after  the  quench  will  grow  exponentially  in  time.  Next,  we 
obtain  a  universal  description  of  interactions  created  by  the 
pseudo-potential  U(r)  by  modifying  the  T-matrix  formalism 
to  take  into  account  Pauli  blocking,  and  apply  these  ideas  to 
the  Stoner  and  BCS  instabilities. 

Cooperon  -  In  this  section,  we  obtain  the  Cooperon,  C,  i.e. 


1/kja 


FIG.  2:  Pairing  instability,  (a)  “Binding  energy”  of  a  Feshbach 
molecule  in  vacuum  and  in  the  presence  of  a  Fermi  sea  (relative 
to  2ef)  as  a  function  of  interaction  strength,  corresponding  to  the 
real  part  of  the  T-matrix  pole  frequency  flq-o  =  Re[wg=o].  Pauli 
blocking  by  the  Fermi  sea  results  in  stronger  binding  across  the  res¬ 
onance.  The  kink  occurs  when  the  pair  becomes  stable,  (b)  Pairing 
rate  and  (c)  rate  of  change  of  Kinetic  Energy  as  a  function  of  interac¬ 
tion  strength  on  the  BEC  side  for  various  temperatures  \T  =  0  (pur¬ 
ple,  solid),  0.12,  0.22,  0.5,  0.66,  0.75  Tf  (red,  dashed)].  Tempera¬ 
ture  is  more  effective  at  suppressing  pair  production  at  larger  values 
of  kf  a  as  the  binding  energy  is  smaller,  thus  the  peaks  in  (b)  and  (c) 
become  sharper  at  higher  temperatures.  The  peaks  in  growth  rate  and 
kinetic  energy  rate  qualitatively  match  experiments  on.  Sharp  onset 
at  kFd  «  1.1  is  expected  to  be  smoothed  by  three-body  processes. 


the  T-matrix  that  takes  into  account  Pauli  blocking  of  states 
by  the  Fermi  sea  (see  Fig.  [3},  which  is  needed  for  an  accurate 
study  of  instabilities  near  a  Feshbach  resonance. 

In  the  center  of  mass  frame  in  vacuum,  the  scattering  of  a 
pair  of  particles  with  identical  masses  to  near  a  wide  Feshbach 
resonance  is  described  by  the  T-matrix  (scattering  amplitude) 

t(E)  =  —  (-  +  i\frrvE\  .  (2) 

47t  \a  ) 

Here,  E  is  the  energy  of  the  scattered  particles  and  the  pseudo¬ 
potential  U(r  —  r')  that  appears  in  Eq.  [T]  is  related  to  the 
T-matrix  via  the  Lippmann-Schwinger  equation.  To  cor¬ 
rectly  renormalize  the  Cooperon,  we  compare  the  Lippmann- 
Schwinger  equations  in  a  Fermi  sea  and  in  vacuum.  For  the 
Cooperon  we  can  not  just  use  the  center  of  mass  frame,  as  the 
Fermi  sea  breaks  translational  invariance.  Therefore,  we  use 
the  laboratory  frame  for  both  to  obtain 

C-^E.  q)  =  r'1  (E  +  2ef  -  q2/4 to) 

f  d3k  nf(f +k)+nf(f-k) 

J(  2tt)3  E  —  ea+k  —  ea_k  ‘  J 

Here,  E  and  q  are  the  center  of  mass  frequency  and  momen¬ 
tum  of  the  pair,  £/  is  the  Fermi  energy,  n F  (k)  is  the  Fermi 
function,  and  ek  =  k2/2TO  —  e/.  Our  approach  is  analogous 


d(KE/eF)/dt 
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to  the  one  used  for  the  Fermi-polaron  problem  Q3- 

Pairing  instability  -  The  Cooperon  enters  the  pairing  sus¬ 
ceptibility  via 


Xpair (q)  =  J  dhdkzGihWq-k^CWx 

x  G(k2)G(q-k2),  (4) 

where  q  stands  for  the  external  frequency  and  momentum 
vector  {E,  q},  dk\  stands  for  dui i  dki/(27r)4,  and  G(k\)  = 
G(w i ,  ki )  is  the  bare  Fermionic  Green  function  in  the  non¬ 
interacting  Fermi  sea  corresponding  to  the  initial  state.  The 
poles  of  the  pairing  susceptibility  correspond  thus  to  the  poles 
of  the  Cooperon,  whose  structure  we  now  explain. 

We  begin  our  analysis  with  the  T-matrix  in  vacuum.  For 
each  value  of  the  scattering  length,  the  T-matrix  has  a  line  of 
poles  on  the  BEC  side  located  at  ojq  =  Qq+iAq  =  —l/ma2+ 
mq2 /4,  corresponding  to  the  binding  energy  of  a  Feshbach 
molecule  with  center  of  mass  momentum  q.  As  a  consequence 
of  energy  and  momentum  conservation  the  pole  frequency  is 
strictly  real  ( Aq  =  0),  indicating  that  a  two-body  process  in 
vacuum  cannot  produce  a  Feshbach  molecule. 

In  the  presence  of  a  Fermi  sea,  the  states  below  the  Fermi 
surface  are  Pauli-blocked,  shifting  the  poles  of  the  Cooperon 
relative  to  the  T-matrix  in  vacuum  in  two  important  ways. 
First,  the  real  part  of  the  pole  Qq,  which  would  correspond 
to  the  binding  energy  of  a  pair  in  the  absence  of  an  imaginary 
part,  uniformly  shifts  down  (see  Fig.  [2^).  This  shift  is  a  result 
of  Pauli  blocking  m,  and  indicates  an  appearance  of  a  paired 
state  on  the  BCS  side  as  well  as  stronger  binding  of  the  pair  on 
the  BEC  side.  Second,  in  the  range  —  oo  <  l/kpa,  <  1.1  the 
pole  acquires  a  positive  imaginary  part  Aq  that  corresponds  to 
the  growth  rate  of  the  pairing  instability.  As  depicted  in  Fig.[T] 
the  growth  rate  of  the  pairing  instability  increases  exponen¬ 
tially  Ag=o  ~  8e^e7r/2fcFa_2  as  one  approaches  the  Fesh¬ 
bach  resonance  from  the  BCS  side,  i.e.  the  growth  rate  of  the 
BCS  pairing  is  equal  to  the  BCS  gap  at  equilibrium  (T8l.  On 
the  BEC  side,  the  growth  rate  continues  to  increase,  reach¬ 
ing  a  maximum  at  kpa  ss  2,  and  finally  decreasing  to  zero  at 
kpa  ~  1.1,  at  which  point  the  Fermi  sea  can  no  longer  absorb 
the  energy  of  the  Feshbach  molecule  in  a  two-body  process. 
Pairing  deeper  in  the  BEC  regime  takes  place  via  the  more 
conventional  three-body  process  and  would  round  the  pairing 
instability  curve  near  kpa  ss  1.1  in  Fig.  |T] 

So  far,  we  have  concentrated  on  the  instability  at  q  =  0. 
We  find  that  q  =  0  is  indeed  the  most  unstable  wavevector 
throughout  the  Feshbach  resonance,  but  the  growth  rate  re¬ 
mains  finite  up  to  q  =  qcut.  Throughout  the  resonance  the 
approximation  qcut  ss  (-\/3/2)(Ag=o works  reason¬ 
ably  well  except  in  the  vicinity  of  kja  ~  2  where  qcat  reaches 
the  maximal  value  for  a  two-body  process  of  2 kj. 

Stoner  instability  -  One  can  expect  that  a  rapid  quench  to 
the  BEC  side  of  the  resonance,  where  interactions  are  strongly 
repulsive,  results  in  an  instability  towards  Stoner  ferromag¬ 
netism.  We  shall  assume  that  right  after  the  quench,  the  atoms 
are  still  in  the  free  Fermi  sea  initial  state  and  the  Stoner  insta- 


FIG.  3:  Diagrams  for  the  vertex  function  T  and  Cooperon  C.  Solid 
lines  represent  bare  fermion  propagators,  dashed  lines  interactions, 
and  wavy  lines  external  sources.  External  legs,  represented  by  gray 
lines,  are  shown  for  clarity. 
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FIG.  4:  Properties  of  growing  collective  modes  in  the  Stoner  insta¬ 
bility  in  3D.  (a)  Growth  rate  of  the  most  unstable  mode  Ag  as  a 
function  of  wavevector  q  for  T  =  0  and  l/kpa  =  0.85  (top  line), 
0.86,  0.87,  ...,  0.93  (bottom  line),  (b)  The  most  unstable  wavevec¬ 
tor  gmax  (blue)  and  the  corresponding  growth  rate  Amax  (red)  vs. 
l/kpa.  A  fit  to  the  mean-field  critical  theory  (y  =  1/2,  z  =  3) 
is  shown  with  black  lines  ED.  (c)  Details  of  the  critical  behavior 
of  gmax  and  Amax  as  a  function  of  distance  from  the  transition  point 
u  =  (l/fc_Fo)c  —  (l/kpa),  (1/ kpa)c  ~  0.94. 


bility  is  competing  with  the  pairing  instability.  Our  goal  is  to 
compute  the  ferromagnetic  susceptibility  using  the  Cooperon 
to  describe  effective  inter-atomic  interactions.  Using  the  full 
Cooperon  allows  us  to  include  three  important  aspects  of  the 
problem:  energy  dependence  of  the  scattering  amplitude  near 
the  Feshbach  resonance;  Pauli  blocking,  which  renormalizes 
the  energy  of  the  virtual  two  particle  bound  states  involved  in 
scattering;  and  Kanamori-like  many-body  screening  0. 

Technically,  we  compute  the  vertex  function  rw  q(wi,  ki), 
which  is  related  to  the  susceptibility  via 

Xfm(o)  =  J  dk\  G(q  +  k\)G(k\)  F^(fci).  (5) 

We  note,  that  the  poles  of  the  susceptibility  and  the  vertex 
function  coincide.  Replacing  the  point  contact  interaction  ver¬ 
tex  by  the  Cooperon  in  an  RPA  type  resummation  of  the  vertex 
function  (see  Fig.  [3] and  Ref.  0)  we  obtain 

I>(fci)  =  1 

+  f  dk2  T q(k2)C (k\  +  k2  +  q)G(k2  +  q)G(k2).  (6) 
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To  compute  the  vertex  function,  a  number  of  approxima¬ 
tions  are  unavoidable.  First,  we  assume  that  q  and  ui  are 
both  small,  which  is  valid  in  the  vicinity  of  the  Stoner  tran¬ 
sition.  Second,  in  the  spirit  of  Fermi  liquid  theory,  we  as¬ 
sume  that  the  most  important  poles  come  from  the  Green  func¬ 
tions,  and  thus  we  replace  G( k2  +  q, tu2  +  w)G( k2,w)  — l 
fFmJ-q.k2<Ha0(5(lk2l  -  M  itS  I  ■  We  then  obtain 
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ili,(k2)C(kl+k2,  w)/qiW(k2),  (7) 


where 


7q,o;(k2)  — 


k%dk2  nF( k2  ^  q/2)  -  np(k2  +  q/2) 


2n2 


£k2-q/2  +  £k2+q/2 


,  (8) 


and  k  indicates  a  vector  on  the  Fermi  surface.  We  thus  make 
the  approximation  that  we  can  replace  ki  and  k2  by  ki  and  k  > 
when  we  evaluate  the  value  of  the  Cooperon.  In  other  words, 
we  assume  that  the  Cooperon  changes  slowly  compared  to 
the  Green  functions.  The  approximation  is  fully  justified  for 
weak  interactions,  where  the  Cooperon  is  momentum  and  fre¬ 
quency  independent,  and  the  vertex  function  matches  the  RPA 
result  EG).  For  strong  interactions,  the  Stoner  instability  is  not 
driven  by  the  pole  of  the  Cooperon,  and  we  therefore  believe 
that  our  approximation  captures  the  essential  physics. 

In  the  range  —0.2  <  l/kFa  <  1.0,  there  is  one  or  more 
lines  of  complex  poles  with  a  positive  imaginary  part  Aq, 
which  corresponds  to  the  Stoner  instability  in  different  chan¬ 
nels  (a  combination  of  momentum  and  orbital  moment).  As 
q  — >  0,  the  different  instabilities  can  be  identified  as  different 
angular  momentum  channels.  Since  magnetization  is  a  con¬ 
served  order  parameter,  in  each  channel  Aq  grows  linearly 
for  small  q.  At  large  q  the  cost  of  bending  the  order  parame¬ 
ter  results  in  the  vanishing  of  Aq  for  q  >  qcut.  In  between, 
Aq  reaches  its  maximum  value  Amax  at  a  wave-vector  gmax 
which  corresponds  to  the  fastest  growing  mode  (see  Fig.  [4]). 

Discussion  -  The  growth  rates  of  the  pairing  instability 
A®£®  and  the  ferromagnetic  instabilities  in  the  various  chan- 
nels  A™x  are  compared  across  the  Feshbach  resonance  in 
Fig.  □  We  have  also  included  the  naive  RPA  estimate  of  the 
growth  rate  of  the  Stoner  instability  in  which  we  have  replaced 
the  Cooperon  by  47t a/m.  From  the  comparison,  we  see  that 
(1)  the  Cooperon  suppresses  the  growth  rate  of  the  ferromag¬ 
netic  instability  but  does  not  eliminate  it,  (2)  the  pairing  and 
ferromagnetic  instabilities  compete  over  a  wide  range  of  in¬ 
teraction  strength  on  both  sides  of  the  resonance,  and  (3)  the 
pairing  instability  is  always  dominant.  Our  results  suggest  that 
even  if  there  is  a  metastable  ferromagnetic  state  ns,  it  prob¬ 
ably  cannot  be  reached  via  dynamic  tuning  of  the  interaction 
starting  from  a  balanced  gas.  However,  for  short  timescales 
~  (A™x)  ~  (A®£q)  ,  both  pairing  and  magnetic  corre¬ 

lations  will  develop  and  may  be  detectable  experimentally. 

Comparison  with  experiment  -  The  maximum  of  the  pair¬ 
ing  instability  in  the  vicinity  of  kFa  ~  2  closely  matches  the 
location  of  the  transition  found  experimentally  ED.  To  inves¬ 


tigate  further,  we  plot  the  pairing  rate  as  a  function  of  inter¬ 
action  strength  for  several  different  temperatures  in  Fig.  |2j). 
The  shape  of  the  pairing  rate  curve,  especially  at  higher  tem¬ 
peratures  looks  qualitatively  similar  to  the  atom  loss  rate  (into 
pairs)  found  experimentally. 

A  fast  rampdown  of  the  magnetic  field  was  used  to  convert 
weakly  bound  molecules  into  strongly  bound  molecules,  and 
the  kinetic  energy  of  the  remaining  atoms  was  measured.  It 
was  found  to  have  a  minimum  at  kFa  ~  2  El-  We  show  that 
this  minimum  can  be  qualitatively  understood  within  our  anal¬ 
ysis  of  the  pairing  instability.  The  energy  of  each  molecule 
produced  is  given  by  ~  Re[wg]  (see  Fig.  [2J1).  The  molecular 
energy  corresponds  to  the  kinetic  energy  of  the  Fermions  re¬ 
moved  from  the  Fermi  sea,  measured  with  respect  to  the  Fermi 
energy.  Thus  the  rate  of  kinetic  energy  change  of  “unpaired” 
atoms  is  ~  (Re[wg]  —  2eF)  x  Im[wq=o]  (see  Fig.  [2J:).  We 
find  that  the  kinetic  energy  minimum  is  in  the  vicinity  of  the 
maximum  of  the  pairing  rate,  in  qualitative  agreement  with 
Ref.  ITTl. 
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